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one with homogeneous and isotropic properties p(t) and another inhomogeneous and anisotropic 
pin(t,r). The rate of expansion of these evolving wormholes is only determined by the isotropic 
and homogeneous matter component p(t). Particularly, we consider a family of exact two- fluid 
evolving wormholes expanding with constant velocity and satisfying the dominant and the strong 
energy conditions in the whole spacetime. In general, for the case of vanishing isotropic fluid pit) 
and cosmological constant A the space expands with constant velocity, and for p(t) — and A / 
the rate of expansion is determined by the cosmological constant. The considered here two- fluid 
evolving wormholes are a generalization of single fluid models discussed in previous works of the 
present authors [Phys. Rev. D 78, 104006 (2008); Phys. Rev. D 79, 024005 (2009)]. 
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I. INTRODUCTION 

Wormhole spacetimes have become one of the most 
popular and intensively studied topics in general relativ- 
ity. Throughout the last decades there has been an accu- 
mulating volume of works on the analytic wormhole ge- 
ometries. The various approaches include both static [l| 
and evolving relativistic versions They principally 
consider static wormhole spacetimes sustained by a sin- 
gle fluid component which requires the violation of the 
null energy condition (NEC), and the interest has been 
focused on traversable wormholes, which have no hori- 
zons, allowing two-way passage through them. These 
hypothetical tunnels in spacetime allow effective super- 
luminal travels, although the speed of light is not locally 
surpassed 0]. 

It is interesting to note that for constructing wormhole 
geometries in general is adopted the reverse approach for 
solving the Einstein field equations. This means that 
one first fixes the form of the spacetime metric (such as 
the redshift and shape functions) and then, by comput- 
ing the field equations, one finds the energy-momentum 
tensor components needed to support such a spacetime 
geometry. The obtained in such a way stress compo- 
nents automatically satisfy local conservation equations, 
by virtue of the Bianchi identities. This reverse method 
helps us to find that a static traversable wormhole vi- 
olates the NEC [3, Q , thus in general relativity an ex- 
otic type of matter is required for sustaining a static 
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traversable wormhole. It is interesting to note that there 
are explicit static wormholes solutions respecting the en- 
ergy conditions in the whole spacetime in Einstein-Gauss- 
Bonnet gravity @. Notice also that in higher dimensions, 
the presence of terms with higher powers in the curva- 
ture provided by certain class of Lovelock theories, allows 
to remove the possibility of violating energy conditions 
even in vacuum, since the whole spacetime is devoid of 
any kind of stress-energy tensor [7[ • 

However, it is well known that in Einstein gravity there 
are nonstatic Lorentzian wormholes which do not require 
WEC violating matter to sustain them. Such worm- 
holes may exist for arbitrarily small or large intervals 
of time |8| . 

On the purely gravitational side, most of the efforts are 
directed to study Lorentzian wormholes sustained by a 
single exotic fluid in classical general relativity. However, 
one can consider also gravitational configurations filled 
with two or more fluids [9(. For example, in cosmology 
such two-fluid models are widely considered today in or- 
der to explain the observed accelerated expansion of the 
Universe (lp| . 

In this paper we shall study evolving wormholes 
sustained by two fluids: one with homogeneous and 
isotropic properties p(t) and another inhomogeneous and 
anisotropic pi n (t,r). The theoretical construction of 
these wormholes will be performed by imposing condi- 
tions on the stress-energy tensor threading the evolving 
wormhole geometry. Specifically, we shall consider the ra- 
dial and tangential pressures of the inhomogeneous and 
anisotropic matter to obey barotropic equations of state 
with constant state parameters. On the other hand, the 
homogeneous and isotropic fluid is taken to be that of a 



2 



perfect fluid described by the energy-momentum tensor 
T a p = (p + p)u a Uj3 - pg a p, (f ) 

where u a is the four-velocity of the fluid, p and p are 
the energy density and the pressure of the cosmic fluid 
respectively 

We shall suppose that the dynamics of the gravita- 
tional fields is governed by Einstein field equations 

R 

Rap — -^g a /3 — —nTa/3 ~ A(? a/ g, 

where k = 8irG and A is the cosmological constant, and 
the evolving wormhole metric will be given by 



ds 2 



z 2 *^dt 2 + a(t) 2 



dr 2 



1 



He) 
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where &(t,r) is the redshift function, a(t) is the scale 
factor of the wormhole universe, b(r) is the shape func- 
tion and dfl 2 = dd 2 + sin 2 9dip 2 . Note that the essential 
characteristics of a wormhole geometry are encoded in 
the spacelike section of the metric ([7]). It is clear that 
this metric becomes a static wormhole if a(t) — > const 
and r) = <&(r) and, as b(r) -> and r) -> it 
becomes a flat Friedmann-Robertson- Walker metric. 

The organization of the paper is as follows: In Sec. 
II we present the dynamical field equations for worm- 
hole models with a matter source composed of an ideal 
isotropic cosmic fluid and an anisotropic and inhomoge- 
neous one. In Sec. Ill some aspects of the geometry of 
the general solution are discussed. In Sec. IV expanding 
wormholes are discussed, and in Sec. V we conclude with 
some remarks. 



II. FIELD EQUATIONS 

Let us now consider the dynamical field equations de- 
scribing evolving wormhole models ([7]). We shall be in- 
terested in studying wormhole scenarios filled with two 
fluids p = p(t) and p in = p in {t,r), where the first cos- 
mic fluid always remains homogeneous and isotropic and 
the other component is in general an inhomogeneous and 
anisotropic fluid. Since we have a spherically symmetric 
space-time, the cosmic fluid p in (t, r) in general may have 
anisotropic pressures, which we shall define as p r (t, r) and 
Pl(t, r) for the radial and lateral components respectively. 
If p r (t,r) = pi(t,r) we have an isotropic inhomogeneous 
pressure. 

Thus, for spherically symmetric spacetimes written in 
comoving coordinates J7J and filled with these two kinds 
of cosmic fluids, the Einstein field equations may be writ- 
ten in the following form: 



3e -2<t>(t,r) H 2 



V 



np in {t,r) + Kp(t) + A, (3) 



-3*(t,r) [ 2 ^ + H 2 } _ + 2e -2^t,r) H ^_ + (4) 



a o\ b — rb' 



r-a- ' dr 

M[2- + H 2 



2e-^) H ^ + __ {2r _ b _ rbr ^ 



2a 2 r 3 
dr 



(5) 




2«-*<-V^f« = 0,(6) 

where it was assumed that the 4-velocity of both fluids 
is the timelike vector u a = (e - ^, 0,0,0), k = 8nG, H = 
a/ a, and an overdot and a prime denote differentiation 
d/dt and dj dr respectively. 

In attempting to find solutions to the field equations we 
first note that Eq. © gives some constraints on relevant 
metric functions which separate the wormhole solutions 
into two branches: one static branch given by the condi- 
tion a — and another non-static branch for d<p/dr = 0. 

In what follows we shall restrict our discussion to non- 
static branch. The condition d<fi/dr = implies that 
the redshift function can only be a function of t, i.e. 
4>(t, r) = f(t) so, without any loss of generality, by rescal- 
ing the time coordinate we can set <f>(t,r) = 0. Thus 
we shall look for the solutions to Einstein equations de- 
scribed by the metric form 



ds 2 



-dt 2 + a{t) 2 



dr 2 



b(r) 



r 2 dn 2 



(7) 



so we have to put $ = into the Eqs. ©-(HI). As 
we shall see below this will imply that the anisotropic 
and inhomogeneous matter component p in (i, r) cannot 
be isotropic. 

Notice that the field equations ©-(JS]) generalize the 
Einstein equations considered in Refs. and [l2j]. 

Thus, in order to quickly solve the field equations we 
shall use the conservation equations = 0. By sup- 

posing that each fluid satisfies the standard conservation 
equation separately we obtain 



dp^ 
dt 



^+3H(p + p)=Q, 
H(3 Pin + Pr + 2 Pl ) = 0, 

2 (P, ~Pr) _ dp r 

r dr ' 



(8) 
(9) 
(10) 



where Eq. ((HJ states the conservation of the isotropic and 
homogeneous component and Eqs. ((9]) and (fT0|) are valid 
for the anisotropic and inhomogeneous cosmic fluid and 
may be interpreted as the conservation equation and the 
rclativistic Eulcr equation (or the hydrostatic equation 
for equilibrium for the anisotropic matter supporting the 
gravitational configuration) respectively. 
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Notice that from equations ([9]) and (fT0|) we see that 
for an isotropic but still inhomogeneous matter com- 
ponent p in , i.e. p x — p r = p in , we have to require 
dp r /dr = 0, so the pressure will depend only on time t, 
obtaining the standard cosmological conservation equa- 
tion p in {t) + 2>H[p in (t) + p i7l (t)} — 0. Thus in this case 
we have a noninteracting superposition of two homoge- 
neous and isotropic fluids p and p in . This leads us to 
conclude that if we want to study evolving wormholes 
filled with a mixture of homogeneous, isotropic fluid p(t) 
and an inhomogeneous fluid p in (t,r), we must consider 
only anisotropic matter component p in (t, r) withp r =/= p l . 

For solving the field equations of the considered gravi- 
tational configuration we shall consider the following an- 
zats: we shall require that the radial and the lateral pres- 
sures have barotropic equations of state. Thus we shall 
write for them 



p r (t,r) = u r p in (t,r), 



(11) 



where w r and lu 1 are constant state parameters (note that 
in this case we have that p g — p^ — p l due to the spherical 
symmetry) . 

In the following, for solving the field equations, we shall 
require that p r (t,r) ^ p ( (i,r). By taking into account 
Eq. ([IT]), from Eq. (dUJ) we get 



p. n (t,r) = F(t) r 2(«,-«r)/«r j 



(12) 



where F(t) is an integration function, and introducing 
this expression into Eq. ^ we have for the energy density 
of the anisotropic matter 



ft„(*,r) 



(13) 



where C is an integration constant. 

Now, by subtracting Eqs. @ and ([5]), and using the full 
energy density (fT3|). we obtain the differential equation 



(14) 



k(cj, - u r ) C r 2 H -<*•)/<*■ 3fo _ r y 



a (3+w r +2 Wi ) 



2a 2 r 3 



It is straightforward to see that in order to have a solu- 
tion for the shape function b = b(r) we must impose the 
constraint 



U) r + 2uj 1 + 1 = 



(15) 



on the state parameters uj r and w, , thus obtaining for the 
shape function 



b(r) = C 3 r 3 - K,Ccu r r- 1/u 



(16) 



where C3 is a new integration constant. Notice that 
the constraint (|15[) implies that the radial and tangen- 
tial pressures are given by 

p r = U r p in , p l = -i (l + Wr)Pi„, (17) 



so the energy density and pressures satisfy the following 
relation: 



P m +p r + ^Pi = 0. 



(18) 



This equation implies that the inhomogeneous and 
anisotropic component satisfies the strong energy con- 
dition. In this case for oj r < —1, —1 < u r < 1 and 
Lu r > 1 we have that uj 1 > 0, — 1 < uj 1 < and uj 1 < — 1 
respectively. 

Now, from Eqs. (j3|), ([T3|) . (fT6|) and taking into account 
the constraint (|15j) we obtain the following master equa- 
tion for the scale factor: 



3H Z = 



3C, 



np - 



A. 



(19) 



Note that by taking into account the metric ([7]) we 
conclude that C 3 may be absorbed by rescaling the r- 
coordinate as follows: C3 = 1 for C3 > and C3 = — 1 for 
C3 < 0, so without any loss of generality we can identify 
it with the spatial curvature parameter k by putting C 3 = 
k, with k = — 1, 0, 1. Thus the master equation (fT9|) ) may 
be written in the form 



3fc 

3# 2 + -5- = np(t) + A. 



(20) 



Summarizing, we have shown that for the gravitational 
configuration 



dr 2 



ds 2 = dt 2 - a(t) 2 x 



+ r 2 (d9 z + sin 2 6d<p 2 ) 



1 — kr 2 + kC UJ r r 



filled with the inhomogeneous cosmic fluid 



Cj .-3-l/u 



(21) 



(22) 



(with anisotropic pressures p r — 0J r p irl and pi — — ^(1 + 
ui r )p i7i ) and another noninteracting arbitrary homoge- 
neous and isotropic p(t), the evolution of the scale factor 
a(t) is governed by the standard Friedmann equation (|20|) 
and the conservation equation ([5]). 

In conclusion, the rate of expansion of these evolving 
wormholes is only determined by the matter component 
p(t) which may be in principle an ideal barotropic fluid, 
a scalar field or any other isotropic and homogeneous 
cosmic fluid considered in literature. Notice that if p(t) = 
A = the space expands with constant velocity 
and for p(t) = and A 7^ the rate of expansion is 
determined by the cosmological constant [l2[ . 

It is worth noticing that for uj r = —1/3 we have a lim- 
iting case since the 3-space becomes isotropic and homo- 
geneous and the anisotropic matter behaves as an ideal 
string gas (i.e. p r = p l = —p in {t)/3). From Eq. (|2"2")l 
we conclude that in this case the matter component be- 
haves as p. n — C/a 2 (t) 7 implying that we have a FRW 
cosmology filled with a mixture of a curvature fluid with 
a cosmic fluid p(t). 
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SOME REMARKS ON THE GEOMETRY 
OF THE SPACE-TIME 



III. 



It is clear that the gravitational configuration ([21]) is 
sustained via a matter source made of the inhomogeneous 
and anisotropic cosmic fluid ([22]) . Let us point out some 
properties of the discussed geometry. In general the met- 
ric (|2ip is not conformally flat since the Weyl tensor does 
not vanish for this metric, except for C = or C ^ and 
u r = —1/3. On the other hand, this acceleration-free 
space-time {g u = 1) is characterized by zero anisotropic 
stress a a p{t,r) and zero heat- flux vector q a . 

Note that the metric ([2~Tj) is conformal to the following 
static metric: 



ds = dr 



dr 2 



1- kr 2 + nCuJrr- 1 - 1 /^ 



+ r 2 dW , (23) 



2 and t = J dt/a(t) is the 



where dfl 2 = d9 2 + sii 

conformal time. In general the component g~^ of the 
metric (|23l) may be valid for all r > or vanish for some 
value of the radial coordinate r > 0; however this does 
not mean that this space-time contains an event horizon 
at ro since g TT = 1. So in principle, for some sets of the 
model parameters, this space-time may contain a naked 
singularity at r = which may be observable from the 
outside. 

In general this geometry admits three-dimensional 
slices t = t a = const with a variable curvature. In this 
case the 3-curvature may be written as 



>R = 



6k + 2kCV 



-3-Tf- 



(24) 



where ao = a(io). If w,. < —1/3 or uj t > these slices are 
asymptotically flat for k = 0, asymptotically de-Sitter 
for k = 1 or asymptotically anti de-Sitter for k = — 1. 
Note that for these ranges of uj r there may arise a naked 
singularity at r — 0, and on the other hand the energy 
density of the inhomogeneous matter component vanishes 
as r — > oo since from Eq. (|2"2"|) we have that p in (to, r) — > 0. 
It is remarkable that the 3-dimensional slices t = to of 
metric (121 1) include as a particular case the 3-dimensional 
slices t = to of the Kottler metric 



ds 2 



1 



2M A 
~r~ ~ 3~ ? 



dt 2 



dr 2 



2M 



At-2 
3 ' 



r 2 dQ 2 , 



(25) 



which includes the de-Sitter (A > 0) and anti de-Sitter 
(A < 0) space-times. The 3-curvature of its slices t = to 
is given by 3 R = — 2A, so they are also spaces of constant 
curvature. By comparing metrics (I21[) and (|25p we con- 
clude that Kottler slices t — to are obtained from slices of 
metric (f2Tj) by putting a,o~ 2 k — A/3, nCu) r = — 2M and 
uj r — > ±oo. From Eq. we see that 



'Rai 



-6k - 2nCr~ 



-6k, 



for ui r — > ±oo, thus effectively we have in this case a space 
of constant curvature. It can be shown that for this limit 
the inhomogeneous energy density ([22]) vanishes since 



Q r -3-l/u 

a 2 {t) 



Curr' 3 - 1 / 1 * 
uj r a 2 (t) 



(26) 



for w r — > ±oo, while the anisotropic pressures take the 



forms p r (t, r) 



'-s and p, (t, r) 



Unfortu- 



nately this model with a vanishing inhomogeneous and 
anisotropic energy density p (n (t, r) and non-vanishing 
pressures p r (t,r) and p,(t,r) is non-physical, so we rule 
it out from consideration. 

If — 1/3 < LJ r < the curvature increases with radius 
and as r — > oo we have the asymptotic forms 3 i? ~ r a , 
p in r-j r a and f 2 (r) ~ with a > and (3 > 2 respec- 
tively. These gravitational configurations do not contain 
any singularity at r — for slices t = to- 



IV. EXPANDING WORMHOLE UNIVERSES 

Now we shall study gravitational configurations where 
the considered solution represents an expanding worm- 
hole geometry [J, |5j . Before treating Lorentzian worm- 
hole geometries in more detail, let us note again that the 
metric ansatz ([7]) provides an explicit class of dynamic 
wormholes that generalize the static, spherically sym- 
metric ones first considered by Morris and Thorne 
Several other aspects of static and evolving wormhole 
spacetimes are analyzed in Refs. [l[ and respectively. 

In order to have a wormhole geometry the functions 
$(r) and b(r) must satisfy some constraints defined by 
the authors of Ref. 4]. However, these authors origi- 
nally were interested only in wormhole geometries fea- 
turing two asymptotically flat regions connected by a 
bridge. In our case we are interested in a more general 
asymptotic behavior. Due to the presence of the constant 
k = —1,0,1; besides the asymptotically flat wormholes 
(k = 0) we may have anti-de Sitter (fc = —1) asymptotic 
wormholes which also may be of particular interest [l3j]. 
Thus the main constraints may be defined as follows: 

Constraint 1: A no-horizon condition, i.e. e*^ 1 *^ is 
finite throughout the space-time in order to ensure the 
absence of horizons and singularities. 

Constraint 2: The shape function b(r) must obey at 
the throat r — ro the following condition: b(ro) = ro, 
being ro the minimum value of the r-coordinate. In other 
words g~r( r o) = 0. 

Constraint 3: Finiteness of the proper radial distance 



l(r) = ± 



dr 



y/1 - b(r)/r 



(27) 



for r > ro throughout the space-time. The ± signs refer 
to the two regions which are connected by the wormhole. 

Let us now consider the possibility of having a worm- 
hole geometry. From the metric ([2~Tj) we conclude that 
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e *(r) _ ^ thus the first constraint is automatically ful- form: 
filled. It must be remarked that, for a general dynamical 
wormhole, the definition of the location of a wormhole 
throat is not straightforward. In this case, the posi- 
tion of the wormhole throat depends on the time-slicing, 
and for a time-dependent wormhole it may not be pos- 
sible to locate the entire throat within one time slice, 
as the dynamic throat is an extended object in space- 
time [llj]. In our case, the evolving metric (|2"Tj) is con- 
formal to the static spacetime (|2"3")l , which represents a 
static wormhole for uj r < — 1 or uj t > 0, and then we al- 
ways may determine the location of the wormhole throat 
on the hypersurface with t = to = const. In general, for 
time-dependent spherically symmetric wormhole space- 
times, alternative definitions of a wormhole throat are 
required, equally valid for static as well as for dynam- 
ical wormholes. Several different definitions have been 
given m Refs. H MMi- 

For example Hochbcrg and 
Visser [l4[ and Hayward [l6[ have introduced two inde- 
pendent quasilocal definitions of a throat for dynamical 
wormholes. These authors do not consider global prop- 
erties of the wormholes, i.e. they make no assumptions 
about symmetries, asymptotic flatness, topology, etc., 
and the wormhole throat is a two-dimensional surface 
of nonvanishing minimal area on a null hypersurface. On 
the other hand, in the Ref. [l5[ the authors have defined 
a wormhole throat quasilocally in terms of a surface of 
nonvanishing minimal area on a spacelike hypersurface. 
Some properties of these three definitions are compared 
in [H. 

(ro) 



From the second constraint, g~ r (ro) = 0, we can find 
the minimum value ro of the radial coordinate where the 
wormhole throat must be located. From this throat con- 
dition, and by taking into account the metric (|2ip . we 
obtain for the integration constant 



_ (fcrg - 1) a +Ur )/u 
o — r 



(28) 



yielding for the shape function b(r) and the metric com- 
ponent g rr 



-1/u 



b(r) — r — 



3 / f N -(l+3w r )M- 



ro J \ V r , 



a (*)ffrr = 1 



r 



'•o 



ro 



-krtl L ) fl-(- 



ro 



-(l+3u! r )/u r 



(29) 



respectively. 

It is easy to verify that the wormhole throat is located 
at ro since 6(ro) = ro. In this case the energy density of 
the matter threading the wormhole takes the following 



1 



r 2 uj r a(t) 2 \r a 



-(l+3w r )/«r 



(30) 



Clearly, in order to have an evolving wormhole we must 
require uj r < -1 or w r > (in both of these cases, 
in the g rr metric component, (1 + u> r )/uj r > and 
(l + 3w, )/o; r > 0), implying that the inhomogeneous and 
anisotropic cosmic fluid (|30j) can support the existence of 
evolving wormholes. It can be shown that the form of the 
wormhole is preserved during all evolution. Notice that 
for an anisotropic matter with u r < —1 we have that 
p in (t, r) > 0, while for an anisotropic matter with co r > 
we have that p in (t, r) < 0. For p(t) = A = the shape 
of the wormhole expands with constant velocity In 
the presence of a cosmological constant with p in (t, r) = 
the wormhole configurations have an accelerated expan- 
sion (contraction) [12( ■ Other properties of such evolving 
wormholes are discussed by authors of the Ref. [TJ El 
and references cited therein. 

As explicit examples of evolving wormholes, let us first 
consider the case where the isotropic and homogeneous 
component is given by a perfect fluid with the barotropic 
state equation p(t) — u>p(t) and k = A = 0. Thus the 
scale factor is given by 



a(t) = Qq t K ' 



and the energy density by 
P(t) = 



3k{1+luY t 2 ' 
while the metric takes the form 



(31) 



(32) 



dr 2 



-(H 



,)/u 



ds 2 = dt 2 - ffl 2 t 4/3(l+ U ) x 



\de 2 



(33) 

In this case, the energy density of the anisotropic matter 
and its pressure components are given by 



np in (t, r) 



wrrgogt 4 /^ 1 " 



ro 



-(l+3uj r )/a» r 



(34) 



and (jTTJ) respectively. Note that if oj r < — 1 or uj t > we 
have asymptotically flat FRW regions for r — > oo and 
u) > — 1. Let us suppose that the isotropic fluid satisfies 
the dominant energy condition, i.e. \p\ < p, p > 0. Thus 
if— 1<cj<— 1/3 the expansion of the evolving worm- 
hole is accelerated, i.e. both universes and the throat 
of the wormhole are simultaneously expanding with ac- 
celeration, while for — 1/3 < uj < 1 the expansion is 
decelerated. 

On the other hand, the total matter content is given 
by p total — p(t) + p in (t, r) and for any oj r < we have 
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that p tota[ > 0. For u r > we can have in general time 
intervals where the total energy is positive or negative. 
For oj > —1/3 the wormhole model starts with a posi- 
tive total energy density (since for a fixed r = const the 
isotropic component dominates over the another one), 
then decreases till zero at certain t eq , and becomes neg- 
ative for t > t eq . For uj < —1/3 the total energy density 
starts negative, then increases till zero at certain t eq , and 
becomes positive for t > t eq . 

The case ui = —1/3 is more interesting, since it allows 
us to consider evolving wormhole models satisfying the 
dominant energy condition (DEC) in the whole space- 
time. In this case the total energy density and the cor- 
responding total pressure components are given by 



p + pi 



(r/roY 



Pi 



u r r 2 a\ 



P r ,T — ~^P + ^rPin, 

1 1 



(35) 



From this expressions we conclude that the total matter 
content satisfies the strong energy condition since p T + 
p rT + 2p lT = 0. 

In order to fulfill the DEC we need to satisfy the fol- 
lowing conditions: p T > 0, p T + p rT > 0, p T —p r T > 0, 
Pt + Pi t — an d P T ~ Pi t — 0- These conditions imply 
that the following constraints must be satisfied: 



1 

3ui r 


^22 
< a r 0: 


(36) 


1 + 0J r 

2co r 


^22 


(37) 


1 — LU r 

4w r 


^22 

< a o r o> 


(38) 


3 + 0J r 

8uj r 


^22 

< a o r o- 


(39) 



Hence, for any given ui r > 0, the DEC is fulfilled by 
choosing the parameters a 2 , and r\ satisfying the condi- 
tions (| 3 6 1) - ([3"9"1> . Note that in this case the isotropic fluid 
gives a scale factor given by a(t) = t. Thus, the class 
of analytic two- fluid evolving wormholes with u> = —1/3 
satisfies the dominant and strong energy conditions in 
the whole spacetime and expands with constant velocity. 
It is interesting to note that the violation or not of the 
NEC, at and near the throat, of a dynamical wormhole, 
may be properly connected to the generalizations of the 
flare-out condition for an arbitrary wormhole discussed 
by above cited authors of the Refs. p^4l6j . In the defi- 
nitions of [Til [iH ] dynamical wormhole throats are trap- 
ping horizons, i.e. hypersurfaces foliated by "marginally 
trapped surfaces", and the violation of the NEC is a 
generic property of such wormhole throats. While in the 
Ref. [151 ] a dynamical spherically symmetric wormhole 
throat is defined as a "trapped sphere", and the NEC 
can still be satisfied for some wormhole configurations. 



As a second example, we shall consider an evolving 
wormhole with k = and filled with a minimally cou- 
pled scalar field with the exponential potential V(<p) — 
Vbe~ A *, where A and Vq are constant parameters. In 
this case the energy density and the pressure of the ho- 
mogeneous and isotropic component are given by p^ = 

2 /2 + V{4>) and p^ = </> 2 / 2 - V(<f>) respectively. Thus 
Eqs. © and imply that 



cj)(t) + 3H<j) - XV o e- A0 = 0, 



(40) 
(41) 



A particular exact solution, describing the power-law ex- 
pansion 



a(t) — aat p 
of the evolving wormhole, is given by 



where 



m = - In t, 



X 2 2k t/ P(3p- 1) 
A = — ,vo = • 

p K 



(42) 



(43) 



(44) 



The first relation implies that p > 0, and from the second 
one we have that Vb<0if0<p<l/3 and Vq > if 
p>l/3. 

Finally, as a last example we shall consider evolving 
wormholes with k = and filled with a tachyon field giv- 
ing the power- law expansion (|42[) , where p is a constant 
parameter. In this case the energy density and the pres- 
sure of the homogeneous and isotropic component are 
v( 



jiven by p T = 



and p T = —V((f>)y 1 — 



tively. Thus Eqs. © and (HOI) imply that 



ZH<t>- 



1 dV(<f>) 
V{$) d<t> 

r2 KV(4>) 



0. 



3ir 



respec- 

(45) 
(46) 



1 



It can be shown that in order to have the power-law 
expansion (|42p the tachyon potential takes the form 
V{4>) — a<j)~ 2 , where <f>{t) = fat, <j)o — \/2/3p and 
net = 2py/l - 2/(3p). 

Notice that for discussed expanding wormholes, filled 
with a scalar and tachyon fields, the wormhole geometry 
is given by 



dr 2 



1- l± 



-(1+Ur)/U, 



ds 2 = dt 2 - a 2 t 2p x 



+ r 2 {d9 2 + sin 2 6dip' z ) 



(47) 
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In this case, the energy density of the anisotropic com- 
ponent and its pressure are given by 



1 



w r r 2 alt 2 v \r 



(48) 



and (p~7|) respectively. 



V. CONCLUSIONS 

We have developed models for evolving wormholes 
sustained by two cosmic fluids: one with homogeneous 
and isotropic properties and another inhomogeneous and 
anisotropic. It is remarkable that the energy density of 
the matter threading and sustaining such a wormhole is 
the inhomogeneous and anisotropic component, while the 
rate of expansion of the evolving wormhole is determined 
by the isotropic and homogeneous component. This mat- 
ter component may be in principle an ideal barotropic 
fluid, a scalar field or any other cosmic fluid satisfying 
the homogeneity and isotropy requirements. For the case 
where the cosmological constant and the isotropic and 
homogenous component are absent the inhomogeneous 
space expands with constant velocity, and when only the 
isotropic and homogenous component is absent the rate 
of expansion is determined by the cosmological constant. 

In general, we have wormhole universes for u> r < — 1 
or w r > 0. The present results generalize our previous 
works [ll[ and [l2j]. The case when the studied inho- 
mogeneous geometry represents wormhole configurations 
expanding with constant velocity, i.e. w r < — I or uj r > 
and p(t) = A = 0, was discussed in Rcf. [ll(, while 
the scenarios where the expansion rate is determined by 
the cosmological constant, i.e. p(t) = 0, were discussed 
in Ref. [12J. If now pit) ^ and A = 0, the expan- 
sion rate of the wormhole is determined by this isotropic 
and homogeneous fluid. It is interesting to note that 
the results of Rcf. were generalized to the case of 
evolving wormholes sustained by a single inhomogeneous 
and anisotropic fluid g(t, r), by imposing the generalized 
equation of state g + aP r + 2f3P t — 0, where a and (3 
are constant parameters, and P r and Pt are the radial 
and transverse pressures, respectively [lj|. All particu- 
lar wormhole solutions discussed in this Ref. are related 
to solutions reported in Ref. fl2| . 

Note that the wormhole geometries (|31[) - (|4"5)) far from 
the throat look like a flat FRW Universe. At first glance, 
if the wormhole throat is located outside of the cosmo- 
logical horizon of any observer, then he is not in causal 
contact with the throat. Thus, for late times, an ob- 
server in this wormhole Universe located too far from 
the wormhole throat will see the Universe isotropic and 
homogeneous and it will be in principle unable for him 
to make a decision about whether he lives in a space of 
constant curvature or in a space of a wormhole spacetime. 

An interesting feature of the discussed here solutions 
is that one can consider evolving wormholes with posi- 
tive total energy density, i.e. p total = p(t) + p in (t, r) > 0. 



We always are free to consider positive homogeneous and 
isotropic energy density p(t). In order to have positive 
p in (t,r) we must choose u> r < — 1. In this case the 
inhomogeneous and anisotropic matter component sus- 
taining the evolving wormhole may be considered a gen- 
eralization of the hypothetical phantom energy used in 
cosmology in order to explain accelerated expansion of 
the Universe. Effectively, this cosmic phantom source 
is characterized by a positive homogeneous energy den- 
sity, i.e. p DE {t) > 0, and by an isotropic pressure sat- 
isfying p DE (t) < -p DB {t). Clearly for p DE = ■yp DE 
we have that 7 < — 1. In our case the phantom en- 
ergy is realized by the inhomogeneous matter component 
p in (t, r) > with linear but highly anisotropic equation 
of state p r < —p in < 0, p, > 0. Spherically symmetric 
distribution of phantom energy, depending only on the 
radial coordinate r, with such linear equation of state for 
the radial pressure were considered in Refs. [f7T |. where 
the authors constructed static wormholes sustained by a 
positive energy density p{r) > 0. 

Lastly, let us note that, to our knowledge, the results 
described in this article, for evolving lorentzian worm- 
holes, with energy-momentum tensor associated with a 
mixture of one isotropic and homogeneous fluid with an 
inhomogeneous and anisotropic component, are firstly re- 
ported here. There are in the literature many works re- 
porting on evolving wormhole solutions to the Einstein 
field equations, however, a large number of these pa- 
pers discuss dynamic wormhole solutions with a single 
fluid source [3, IS, uM, LLSj ■ It is worth to mention here 
that a general class of higher evolving dimensional worm- 
holes sustained by a single fluid was studied in Ref. [2(| ■ 
Most specifically, it was considered a quasi-static spher- 
ically symmetric evolving wormholes, with static four 
non-compact dimensions and an arbitrary number of ex- 
tra time-dependent compact dimensions. The results of 
the study show that the WEC cannot be satisfied at the 
throat. This is mainly due to that the matter content is 
not distributed in the whole space and is matched to the 
vacuum. The presence of this matter-vacuum boundary 
places restrictions on the time dependence of extra com- 
pact dimensions, consequently implying the violation of 
the WEC. On the contrary, as shown in the previous Sec, 
we can have evolving wormhole configurations satisfying 
the WEC. 

On the other hand, with respect to wormholes involv- 
ing a mixture of two fluids, a dynamical wormhole, filled 
with a perfect fluid and a ghost scalar field, is provided in 
Ref. [la ]. The considered in this Ref. dynamical worm- 
hole metric (4.53), 



ds 2 = dt 2 



a 2 {t) (dx 2 - (x 2 + b 2 )dn 2 ^J , 



may be rewritten, by using the transformation x 2 = r 2 — 
b 2 , as 



ds 2 =dt 2 -a 2 {t) 



dr 2 



l-b 2 /r 2 



r 2 dn 2 
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where b is a constant parameter and a (t) = t/to- Thus 
the solution (4.58)-(4.60) of the Ref. [l5| is a particu- 
lar case of the discussed in this paper wormhole geome- 
tries (l3"T1) - (j3l?1) . with ui = —1/3, ui r = 1, ro = b and 
Oq = l/io< I n this case the anisotropic and inhomo- 
geneous component with ui r = 1 may be identified as 
a massless ghost scalar field (note that the energy den- 
sity (|34p becomes negative). The conditions (|36l) - (j3T))) 
imply the inequalities Oq7q > 1/3, a^r^ > 1, clqTq > 
and agTg > 1/2 respectively. Hence, for this particu- 
lar solution the DEC is satisfied in the whole space for 
do ro > 1. Notice that the wormhole throat definitions of 



Hochberg-Visser [14| or Hayward [16| do not apply to this 
evolving wormhole since the whole spacetime is foliated 
by trapped surfaces and there is no trapping horizon fljl ] . 

As far as we know, there is only one more paper where 
two-fluid evolving lorentzian wormholes are considered. 



In Ref. [21| FRW models with a traversable wormhole 
are considered. In this case the matter content is di- 
vided into two parts: the cosmic part p that depends on 
cosmological time only, and the wormhole part p w that 
depends on the radial coordinate only. This wormholes 
can be finally connected with the particular wormholes 
solutions (l3T1)-(l34l). 
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